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Influence of the phonon temperature 
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The magnetization response of bulk ferromagnets after excitation by an ultrashort optical pulse is calculated 
using a dynamical model of the Elliott- Yafet type that includes the effects of the spin-orbit interaction in the 
ab-initio ferromagnetic band structure, the electron-phonon interaction at the level of Boltzmann scattering 
integrals, and dynamical changes in the temperature of the phonon bath. Using realistic parameters for the 
ultrashort optical pulse, the computed maximum magnetization quenching achievable with electron-phonon 
scattering in a fixed band structure is much smaller than the quenching observed in experiments. Heating of 
the phonon bath is found to not appreciably change the magnetization dynamics on ultrashort timescales. 
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I. INTRODUCTION 

The magnetization dynamics of ferromagnets after ul- 
trashort pulse excitation has been investigated experi- 
mentally and theoretically since its first observation 15 
years ago.'^ Although a lot of progress has been made 
and a viable macroscopic model for the magnetization 
dynamics exists in the three-temperature model, there 
are still open questions regarding the underlying micro- 
scopic mechanism of the demagnetization as well as the 
details of the angular momentum transfei^^^ on ultrashort 
timescales. Additional interest in this topic is generated 
by its potential implications for laser-induced control of 
the magnetization and, possibly, all-optical switching. 
Some promising microscopic explanations of magnetiza- 
tion dynamics involve mechanisms of the Elliott- Yafet 
type,'^ in which the total spin polarization is reduced by 
scattering processes that change the average spin of the 
scattered electrons. In the presence of spin-orbit interac- 
tion this may be achieved by electron-phonon/impurity 
scattering,^ 5 jg usually assumed, but also by electron- 
electron,!^ electron-magnon,'^ scattering. 

To clarify the importance of electron-phonon scatter- 
ing to the ultrafast demagnetization of ferromagnets, we 
recently investigated this scattering mechanism in bulk 
ferromagnets. By solving the dynamical equations for 
the momentum-resolved carrier distributions we showed 
that the electron-phonon scattering dynamics in the fixed 
band structure cannot explain the observed magnetiza- 
tion dynamics for realistic pulse intensities.'^ Apart from 
the fixed band structure, our earlier calculation used a 
bath assumption for the phonon system and kept the 
lattice temperature fixed at the value before the exci- 
tation. After excitation with appreciable laser intensi- 
ties, however, the electrons should transfer energy to the 
lattice quickly, leading to a rise in lattice temperature. 



The heating of the lattice should, in turn, affect the elec- 
tronic dynamics because the available phase space for 
transitions in the phonon system is changed. Here, we 
investigate this effect on the electronic dynamics by in- 
cluding the temperature of the phonon bath as an addi- 
tional dynamical variable. As in Ref. |8l the optical exci- 
tation and the electron-phonon scattering are described 
at the level of Boltzmann scattering integrals, and the 
band structure, electron-phonon interaction matrix ele- 
ments, as well as dipole matrix elements are calculated 
from density-functional theory. 

II. MODEL 

A detailed account of our numerical method including 
the connection to ab-initio calculations can be found in 
Ref. m In the following, we give a brief description of our 
theoretical approach, and show how a dynamical change 
of the phonon bath temperature is achieved in the frame- 
work of our model. We describe the ferromagnetic dy- 
namics by momentum-resolved carrier distributions n'i (t) 
for all bands fi of interest, for which we integrate the dy- 
namical equations 

dt k dt fcl°Pt dt fcl<=-P ^ ^ 

with the following contributions from the optical excita- 
tion and electron-phonon scattering: 
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Here, E{t) denotes the amplitude of the laser field. 
The scattering kernels /C°p' and /C°"P contain all the 
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occupation independent quantities for the optical and 
the scattering contribution, respectivelyP The scatter- 
ing kernels /C®"P include the electron-phonon matrix ele- 
ments, phonon and electron energy dispersions, and the 
phonon distributions. If one employs a bath approx- 
imation for the phonons and keeps the phonon tem- 
perature Tp fixed, these scattering kernels are indepen- 
dent of time. In this paper, we still assume that the 
phonon system can be treated as a bath, but we char- 
acterize it by a time-dependent temperature Tp = Tp{t). 
This leads to an implicit time dependence of the scat- 
tering kernel /C'P in Eq. ([s]) because it depends on the 
phonon temperature via the phonon occupation numbers 
= l/[cxp{huj~/kBTp) — 1]. Here, ks is Boltzmann's 
constant, and the phonon modes are labeled by their 
polarization A and wave vector q. The computation of 
/C°"P(rp) for a given temperature of the phonon bath im- 
plies a sum over all allowed transitions weighted by the 
corresponding phonon occupation numbers. Instead of 
redoing this summation in Eq. ([s]) in every numerical 
time step, which would vastly increase the CPU time de- 
mands, we note the smooth, monotonic dependence of 
on Tp, which suggests also a smooth dependence of their 
weighted sum /C''"P(Tp). We therefore approximately re- 
place IC°~P{Tp) by a polynomial of degree P which can 
be obtained from a least squares fit to values JC°~^{Tp) 
precomputed for a "grid" of temperatures in the temper- 
ature interval of interest /C'=-P(rp) w YZi^q^TIT'^ ■ 

The dynamical equation of the temperature Tp of the 
phonon bath is derived from the requirement of energy 
conservation, which holds for the system of electrons and 
phonons if one neglects losses to the environment. That 
is, a change in energy of the electronic system E^, should 
lead to an equal, but opposite, change in the energy of 
the phonon system i?p. 
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where is the energy of single-particle state \fi, k) and 
{d / dt)n'i\e-p the scattering contribution to the occupa- 
tion change as in Eq. (|3|. Finally, the phonon heat ca- 
pacity 
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is computed using the phonon frequencies obtained 
from an ab-initio calculation.El As Cp{Tp) also smoothly 
depends on the temperature Tp, it is replaced in the nu- 
merical calculations by an polynomial approximation for 
the same set of "grid" points as /C°"P(Tp). The combined 
set of dynamical equations for the electronic distribu- 
tions, Eq. ([T]), and the temperature Tp of the phonon 



bath, Eq Q, is solved numerically. From the solution, 
the magnetization of the system can be obtained by the 
sum 
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over the spin expectation values (S) ^ j: of the states fc), 
neglecting orbital contributions to the magnetization. 
Here fiB denotes the Bohr magneton. 



III. RESULTS 

We apply the model and calculation procedure de- 
scribed in the previous section to the simple ferromagnets 
nickel and iron, and compare with our earlier result 
that assumed a constant temperature of the phonon 
bath. The parameters were chosen as in Ref. jB) For 
the polynomial interpolations, we evaluated IC°~^ and Cp 
for 40 equidistant temperatures in the range from 250 to 
1000 K, and used an interpolation with third and fifth 
degree polynomials, respectively. 

The results shown in Figs. 1 and 2 for nickel and iron, 
respectively, are computed assuming the same pulse in- 
tensity of 4 mJ/cm^, from which the amplitude of the 
electric field in the material is obtained by material spe- 
cific values for the complex refractive index. The op- 
tical excitation leads to a deposition of energy in the 
material that is quite similar for both nickel and iron, 
cf. Figs. 1(b) and 2(b). The experimentally observed dif- 
ferences between nickeP and irorP at the same pump flu- 
ence are therefore not simply due to a different "net" en- 
ergy transfer to the ferromagnet by the ultrashort pulse. 
The magnetization dynamics exhibits a rapid drop with 
different time constants for iron and nickel of less than 
250 fs, and a subsequent return to an equilibrium state on 
a time scale of 10 ps. The demagnetization occurs via the 
so-called the Elliott- Yafet mechanism: The spin mixing 
due to the spin-orbit interaction leads to a change of the 
carriers' average spin in each scattering event and thus 
to a change of magnetization. On the longer time-scale 
the carrier distributions approach again an equilibrium 
Fermi-Dirac form and thus lead to a return of the mag- 
netization towards an equilibrium value. 

As stressed in Ref. El the computed magnetiza- 
tion dynamics should be compared to experimental re- 
sult^^'^ on a sub-picosecond timescale, because changes 
in the quasi-equilibrium magnetization (exchange split- 
ting), which dominate the dynamics on a picosecond 
timescale, see, e.g., Ref . O are excluded from our model 
due to the fixed band structure. We also neglect carrier- 
carrier scattering which leads to an additional relax- 
ation towards quasi-equilibrium distributions. The com- 
bination of "magnetic" dynamics and electron-electron 
scattering likely prevents the occurrence of a dynamical 
magnetization larger than its equilibrium value, as it is 
the case for nickel around 4ps in our simulation. Figs. 
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1(a) and 2(a) indicate that the dynamical changes in the 
temperature of the phonon bath have no pronounced in- 
fluence on the ultrafast magnetization dynamics. Still, 
the calculation including a time-dependent temperature 
of the phonon bath yields a slightly accelerated magne- 
tization dynamics on the ultrafast timescale because the 
heating of the phonon bath increases the scattering rates. 

Our calculated results on a timescale of a few picosec- 
onds should be interpreted strictly as a model study. 
Figs. 1(b), (c) and 2(b), (c) show that the "heating dy- 
namics" of nickel and iron are very similar, even though 
the magnetization dynamics are distinctly different. In 
both materials, the system reaches a steady state on a 
timescale of several (ten) picoseconds when the tempera- 
ture of the phonon bath is allowed to change dynamically. 
This steady state is characterized by a heated phonon 
system and an higher energy (density) of the electron 
system, because system does not reach the initial equi- 
librium distribution at Tq = BOOK, but instead one at the 
higher final temperature, which corresponds to a higher 
energy (density). For nickel, as shown in Fig. 1, the dy- 
namical phonon temperature even rises above the Curie 
temperature (Tc(Ni) = 631 K)p2!, xhe deposited energy 
of about 130meV/cell is therefore enough to completely 
demagnetize the material after thermalization. Such a 
full demagnetization is not observed in experiment at the 
given intensity, which indicates that the constant ampli- 
tude of the optical field assumed by us is an overestima- 
tion. In the experiment the field is attenuated due to 
absorption in the sample. 



IV. CONCLUSIONS 



ing of ultrafast demagnetization dynamics in ferromag- 
nets. 




Figure 1. (a) Magnetization dynamics after optical excita- 
tion for nickel computed with a dynamical (solid) and a fixed 
(dashed) temperature of the phonon bath; (b) energy in the 
electronic system, and (c) temperature of the phonon system. 



We presented an ab-initio based calculation of electron- 
phonon scattering dynamics with a dynamical tempera- 
ture of the phonon bath for nickel and iron, which con- 
stitutes a more realistic description of a thin film weakly 
coupled to a substrate than the fixed temperature of the 
phonon bath used in a previous study. ^ In the frame- 
work of this model, which uses a fixed band structure, 
we compared the resulting demagnetization scenario with 
constant and time-dependent temperature of the phonon 
bath to estimate the infiuence of the phonon temperature 
on the dynamics. We find only marginal changes in the 
magnetization dynamics. This supports our earlier con- 
clusion- that electron-phonon scattering in a constant 
band structure is not able to explain the pronounced ob- 
served demagnetization for realistic laser intensities. Our 
results provide a hint that a dynamical phonon-bath tem- 
perature is unlikely to play an important role in the ul- 
trafast magnetization dynamics of ferromagnets. Based 
on our results, we conclude that a dynamical change of 
the band structure, i.e., the exchange splitting has to be 
taken into account to achieve a quantitative understand- 



REFERENCES 



^E. Beaurepaire, J.-C. Merle, A. Daunois, and J.-Y. Bigot, Phys. 

Rev. Lett. 76, 4250 (1996). 
2C. Stamm, T. Kachel, N. Pontius, R. Mitzner, T. Quast, K. 

Holldacli, S. Khan, C. Lupulescu, E. F. Aziz, M. Wietstruk, H. 

A. Diirr, and W. Eberhardt, Nat. Mater. 6, 740 (2007). 
3d. Steiauf, C. Illg, and M. Fahnlc, J. Mag. Mag. Mat. 322, L5 

(2010). 

^B. Koopmans, G. Malinowski, F. Dalla Longa, D. Steiauf, M. 

Fiihnle, T. Roth, M. Cinchetti, and M. Aeschlimann, Nat. 

Mater 9, 259 (2010). 
^D. Steiauf and M. Fahnle, Phys. Rev. B 79, 140401(R) (2009). 
"M. Kraul3, T. Roth, S. Alebrand, D. Steil, M. Cinchetti, M. 

Aeschhmann, and H. C. Schneider., Phys. Rev. B 80, 180407 

(2008). 

^E. Carpene, E. Mancini, C. Dallera, M. Brenna, E. Puppin, and 

S. De Silvestri, Phys. Rev. B 78, 174422 (2008). 
*S. Essert and H. C. Schneider, cond-mat:1108.4454 (2011). 
^M. Djordjevic, M. Liittich, P. Moschkau, P. Guderian, T. 

Kampfrath, R. G. Ulbrich, M. Miinzenberg, W. Felsch, and J. 

S. Moodera, phys. stat. sol (c) 3, 1347 (2006). 
"D. E. Gray, ed., American Institute of Physics Handbook, 3rd ed. 

(McGraw-Hill, 1972). 



1.01 n — I — I — 1 — I — I — I — I — r 




2 4 6 8 10 

time (ps) 



Figure 2. Magnetization dynamics after optical excitation for 
iron computed with a dynamical (solid) and a fixed (dashed) 
temperature of the phonon bath; (b) energy in the electronic 
system, and (c) temperature of the phonon system. 



